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Abstract 

In minimal Supergravity (mSUGRA) models the lightest supersymmetric particle (assumed 
to be the lightest neutralino Xi) provides an excellent cold dark matter (CDM) candidate. 
The supersymmetric parameter space is significantly reduced, if the limits on the CDM relic 
density QcDMh 2 , obtained from WMAP data, are used. Assuming a vanishing trilinear 
scalar coupling Aq and fixed values of tan/3, these limits result in narrow lines of allowed 
regions in the m — ui\i% plane, the so called WMAP strips. In this analysis the trilinear 
coupling A has been varied within ±4 TeV. A fixed non vanishing A value leads to a shift 
of the WMAP strips in the m — m^j plane. 



1 Introduction 



The supersymmetric (SUSY a ) parameter space in minimal Supergravity (mSUGRA) scenar- 
ios is usually studied in terms of the common scalar mass mo, the common gaugino mass 
mi/2, the ratio of the Higgs expectation values tan/3 and the sign of the Higgsino mass pa- 
rameter /i. However, the fifth free parameter, the common trilinear scalar coupling A , was 
usually set to zero. In recent studies, the impact of A on the SUSY parameter space was 
recognised [1]. 

In the mSUGRA framework the lightest neutralino lends itself as an excellent cold dark 
matter (CDM) candidate, thus providing a connection between particle physics and astro- 
physics. The inclusion of cosmological experimental data allows to significantly reduce the 
mSUGRA parameter space. The satellite born detector WMAP b measured the abundance 
of CDM in the Universe to be 0.094 < n CDM h 2 < 0.129 (at 2a C.L.) [2]. Under the as- 
sumption of A = TeV and fixed tan/3, only some narrow lines in the mo — mi/2 plane are 
left over as allowed regions after including WMAP data [3]. In this analysis we have studied 
the effects of non vanishing couplings Aq systematically and found extended areas in the 
mSUGRA parameter space, which no longer can be excluded. In the next section we give a 
short theoretical introduction to SUSY and mSUGRA with particular emphasis on the role 
of the trilinear scalar coupling Aq. In section El we describe the constraints from cosmological 
data on the mSUGRA parameter space. The tools used for this work are briefly discussed in 
section El including comparisons of the different Monte Carlo (MC) programs. In section 
we present a scan over the mSUGRA parameter space. Allowing A$ to vary, the WMAP 
strips are broadened to extended areas. A general parametrisation of these allowed regions is 
complicated. Thus, we fixed A to several different values and constructed parametrisations 
for the resulting lines, which are presented in section |H1 

2 Super symmetry 

The Standard Model of particle physics (SM) is in stupendous agreement with experimental 

measurements. Why should it then be extended? The SM encounters several theoretical 

problems, which cannot be solved without the introduction of new physics, (i) In the SM 

the electroweak (EW) symmetry has to be broken in order to generate the masses of the 

weak gauge bosons. The Higgs sector has been introduced to mediate EW spontaneous 

a In the following the shortcut SUSY is used for supersymmetric as well as Supersymmetry. 
k http : // map . gsf c . nasa . gov 
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symmetry breaking. However, the Higgs mechanism is not experimentally established yet. 
(ii) The couplings of the three gauge interactions do not unify at some high energy scale, 
so that the SM cannot easily be included in a grand unified theory (GUT), (in) Due to 
quadratically divergent contributions to the Higgs boson mass, the huge gap between the 
EW and the GUT scale requires the introduction of a fine tuned mass counter term in order 
to establish an intermediate Higgs mass. This problem is known as the hierarchy problem [4]. 
(iv) Furthermore the SM does not provide a candidate for CDM. 

SUSY [5, 6] is one of the best motivated candidates for physics beyond the SM. It cannot 
solve all the problems of the SM but all those listed before. New particles at the TeV scale 
modify the /^-functions of the three gauge couplings such that the latter meet at about 
10 16 GeV [7]. SUSY-GUTs generate the EW symmetry breaking dynamically, if the top 
mass ranges between about 100 and 200 GeV [8], in agreement with the measured top mass 
of 172.7 ± 2.9 GeV [9]. By connecting fermions with bosons the hierarchy problem is solved. 
The quadratic divergences are cancelled systematically order by order, if the corresponding 
couplings between SM and SUSY particles are identical [4] . Fine tuning of the counter terms 
is not required, if the masses of the SUSY particles are not too large, i.e. of the 0(TeV). The 
lightest supersymmetric particle (LSP), if stable, provides a good candidate for CDM [10]. 
However, several problems of the SM remain unexplained in SUSY extensions of the SM as 
e.g. the masses of the fermions or the origin of three generations. 

2.1 MSSM 

The minimal supersymmetric extension of the SM (MSSM) [6], requires a doubling of the 
SM degrees of freedom (d.o.f.) including two complex Higgs doublets H d = (H®,Hj) and 
H u = (H+,H®) giving mass to the down- type and up- type fermions, respectively. For each 
fermionic d.o.f. a corresponding bosonic one exists and vice versa, as can be inferred from 
Tab^ Two complex Higgs doublets are needed for the theory to be supersymmetric and 
free of anomalies. After spontaneous symmetry breaking five of the eight states remain as 
physical particles: two neutral CP-even (scalar), one neutral CP-odd (pseudoscalar) and two 
charged Higgs bosons. The scalar superpartners of the left /right-handed fermion components 
Ql,u r ,cI r and L L ,e R mix with each other yielding the mass eigenstates Qi,2 and L 12 , 
respectively [6]. Since the mixing angles are proportional to the masses of the ordinary 
fermions, mixing effects are only important for the third-generation sfermions. The four 
neutralinos Xi 4 are linear combinations of the SUSY partners of the neutral gauge bosons 
W°, B° and the neutral Higgsinos H° d , the superpartners of the neutral components of the 
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Gauge Multiplets 


Chiral Multi] 


Diets 


J=l 


3=1/2 


3=1/2 


J=0 


Gluon g 


Gluino g 


Quark Q 


Squark Qi y2 


W ± } W° 


Wino W ± , W° 


Lepton L 


Slepton Lx t 2 


B° 


Bino B° 




Higgs H d , H u 



Table 1: Particle content of the MSSM. The gauge multiplets mediate the interactions while 
the chiral multiplets contain the matter content. J denotes the spin quantum number. 

two Higgs doublets. Analogously the charged winos W ± and the charged Higgsinos H£, 
build up the four charginos xt 2 as mass eigenstates. 

The R-parity [11], defined as R = (_ i) 3B + 2S, + L with B = baryon number, S = spin and 
L = lepton number is introduced as a new discrete symmetry, which distinguishes SM 
particles (R = 1) and their SUSY partners (R = —1). In R-parity conserving models the 
SUSY particles (sparticles) can only be produced/annihilated in pairs, so that the LSP is 
stable. 

2.2 mSUGRA 

Since no SUSY particle with the same mass as its SM partner has been discovered, SUSY 
has to be broken. Different scenarios for SUSY breaking mechanisms have been proposed. 
It is typically assumed, that the breaking takes place at a high energy scale. There are 
several models with different messenger particles (gravitons, gauge bosons, . . . ) mediating 
the SUSY breaking effects down to the EW scale. We concentrate on mSUGRA models with 
the graviton as the messenger particle [12]. The key point of these models is the unification 
of the bosonic masses to the common scalar mass mo, of the fermionic masses to the common 
gaugino mass m 1( / 2 and of the trilinear scalar couplings to A at the GUT scale in addition 
to gauge coupling unification. As a consequence, the whole MSSM can be described by only 
five additional parameters to the SM ones (instead of more than 100 parameters as in the 
general MSSM): m , m^ 2) tan/3 and sign(//), the sign of the Higgsino mass parameter. 
All parameters of the MSSM can be derived by renormalisation group (RG) equations from 
the values of these five input parameters (which define a specific mSUGRA model) at the 
GUT scale, as illustrated in FigHJ 

In R-parity conserving mSUGRA models the LSP is a neutralino in large regions of the 
parameter space, if cosmological bounds are included [2]. As the Xi is electrically neutral, it 
does not directly couple to photons, an essential condition for any dark matter candidate. 



3 



4(1(1 



300 



100 

















tan (3 = 1 
Y b = Y T 


65 




















h. \\ 


\uiuino 




m i 










































wino 








in 


n 


1, _Ju 














i 






Mint 












i , , , 


, , , 


, , , 


, , , 




, , , 


, , , 





Figure 1: Unification of the sparticle masses at the GUT scale [13]. The common gaugino 
mass mi/2 and the common scalar mass itlq are input parameters of mSUGRA. 

2.3 Trilinear scalar coupling Aq 

In Fig 121 the linear A dependence of the trilinear couplings of the third generation, A tt t, jT , 
is depicted. The evolution of A down to the EW scale determines the couplings Ai,<2,z f° r 
up-type, down-type squarks and charged sleptons resulting in the approximate relation: 

A k = d k A + d' k my2 with k = u,d, I. (1) 

The coefficients d k depend on the corresponding Yukawa couplings and are of 0(1), increasing 
for decreasing masses. As the top quark is much heavier than the bottom quark and the 
tau lepton, the slope of A t is smaller for small tan/3 (FigEJ. The coefficients d' k depend 
additionally on the gauge couplings and are of order unity [15]. The RG evolution of A Utd j 
from the GUT to the EW scale generates the corresponding terms of the SUSY breaking 
part of the Lagrangian in the low energy limit: 



soft fc\-\ IT % 3 



V2Mi 



w 



%A u HlQi L u\ + ^-A d H\Qi L d\ + ^AH d L{e\ 



sin/? " « 1 cos/3 " _l h CQs/3 



(2) 



where g denotes the coupling constant of SU(2) l and is the two-dimensional antisymmetric 
tensor. M w is the mass of the W boson and M u ^i are the masses of the up/down- type quarks 
and the charged leptons, respectively. These terms of the Lagrangian introduce interactions 
between MSSM Higgs bosons, "left-" and "right-handed" sfermions, in addition to the usual 
interaction mediated by the Yukawa couplings. These additional couplings are proportional 
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Figure 2: Dependence of the trilinear scalar couplings A t , Aj, and A T at the EW scale on A 
for the Snowmass point SPS1A [14]. Values of A outside the displayed region are excluded 
because of numerical problems of the RG equations. 

to the masses of the corresponding fermions (EqJ2J), so that they are only relevant for the 
third generation. 



3 Relic density 

In the early Universe the interaction rate (T) of a species of particles must be larger than the 
expansion rate of the Universe. Otherwise the thermal equilibrium cannot be maintained 
and the particles decouple at the freeze out temperature Tf, i.e. if the following condition is 
fulfilled: 

r = n (cry) = H at Tf, (3) 

where n denotes the number density of the particle species, H is the Hubble constant and 
(ay) is the thermal average of the total annihilation cross section times the velocity of the 
corresponding particles. The ratio of the number of relic particles to the total entropy in 
the Universe remains constant after they are frozen out. Since the present entropy density 
is known, the present number of these particles can be approximately determined by using 
the freeze out condition (Eq|21). The relic density of a particle species X can be estimated 
as [16]: 

„ , 3x KT 27 cm 3 ^ 1 

tt x h 2 « j— , 4 

{o-vjx 
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Figure 3: Difference between the values of the relic density for calculations without (blue 
circles) and with (red points) coannihilation processes. The input parameters are taken from 
the Snowmass point SPS1B (m = 200 GeV, my 2 = 400 GeV, tan/3 = 30, fi > 0) and 
A has been varied within ±1 TeV. The green shaded area shows the region allowed by the 
WMAP data. 



where h denotes the Hubble constant in units of 100 km sec" 1 Mpc" 1 . The cross section 
is proportional to the squared matrix element. The matrix element for a given process is 
inversely proportional to the squared mass of a heavy propagator particle. Assuming SUSY 
CDM, the propagator particle may also be a sparticle. In the case of neutralino annihilation, 
the cross section contains several different channels, thus the dependence of Q x h 2 on the 
mSUGRA parameters is not trivial. 

Moreover, it is important to include all possible coannihilation processes [17] between the 
LSP and the next heavier sparticles as can be inferred from FigOlfor the Snowmass point 
SPS1B (mo = 200 GeV, my 2 = 400 GeV, tan/3 = 30, Li > 0) [14]. A has been varied within 
±1 TeV c . Without coannihilations the Snowmass point SPS1B would be excluded by the 
WMAP data. Including coannihilations several models with A Q ^ TeV are allowed. The 
formulas including coannihilations are more involved [17]. 

The relic neutralino density (incl. coannihilation) for tan (3 = 10, Aq = — 500 GeV and 
c Numerical problems of the RG equations occur for A < — 350 GeV and A > 720 GeV. 
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Figure 4: The relic density Q x h 2 as a function of my 2 obtained using the ISAJET and the 
DarkS US Y Monte Carlo programs, for tan/3 = 10, fi > 0, my 2 < 1 TeV, mo varied between 
and 300 GeV and A = -500 GeV (left) and A = 2000 GeV (right). The red lines belong 
to mo = 190 GeV (left) and = 235 GeV (right). Each single point denotes one set of 
mSUGRA parameters. The green shaded area shows the region allowed by the WMAP data 
and the yellow shaded area indicates the resulting allowed my 2 region. 

+2000 GeV as a function of my 2 with < m < 300 GeV is depicted in FigEJ The lines 
result from fixed values of mo, the step sizes for mo,™^ are chosen to be 5 GeV for both 
plots. This is exemplified by the red lines, which correspond to mo = 190 GeV in the left 
plot and to mo = 235 GeV in the right plot in Fig0] Values below the red lines belong to 
smaller mo values and lines above to larger ones. Models with my 2 170 GeV in the left 
plot and with myi < 480 GeV in the right plot are excluded by collider constraints. With 
growing my 2 the lower limit of the relic density becomes larger. 

Assuming that CDM entirely consists of LSPs, the application of cosmological constraints 
on the mSUGRA parameter space is possible. The WMAP experiment constrained the relic 
CDM density to the narrow range 0.094 < QcDAih 2 < 0.129. For vanishing Aq, only some 
narrow strips remain as can be inferred from FigJSJ The left plot has been obtained by the 
authors of Ref. [3] using their own Monte Carlo programs. The parametrisation of these data 
defines the WMAP strips (black lines in the right plot). The right plot displays the results 
of this study, using the ISAJET and DarkSUSY Monte Carlo generators, for comparison. 
For small values of tan f3 the agreement is good, but the differences increase with tan (3. The 
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Figure 5: Regions in the tuq — my 2 plane accounting for 0.094 < ^\^ 2 < 0.129 for different 
tan/3, n > and A Q = 0. Left plot: WMAP strips from Ref. [3] are shown with the post 
WMAP benchmark points (A \ . . M'). Right plot: The corresponding allowed regions obtained 
with ISAJET and DarkSUSY are shown. The black lines are the parametrisation given for 
the WMAP strips in the left plot. 



mSUGRA models with large values of tan (3, m and my 2 , shown in the left plot (in vicinity 
of M'), are excluded in the right plot since the relic densities are too large. The extended 
region for tan/3 = 50 and small values of m , my 2 obtained with ISAJET does not show 
up in the left plot. If we include variation of the parameter A within ±4 TeV, the WMAP 
strips are broadened to extended will be shown in section El 



4 Monte Carlo generators 

As already stated in section El the RG equations have to be applied to evolve the mSUGRA 
parameters from the GUT scale down to the EW scale. For these evolutions two different 
Monte Carlo programs have been used: 

1. SuSpect [18]: SuSpect 2.2 is a Fortran code used to determine the SUSY particle 
spectra. The calculations are done in the MSSM framework assuming R-parity and 
CP conservation. It can also be applied to constrained scenarios as mSUGRA, anomaly 
mediated SUSY breaking (AMSB) and gauge mediated SUSY breaking (GMSB) mod- 
els. The algorithm incorporates RG equations to evolve the parameters between the 



8 



EW and the GUT scale complemented by constraints from radiative EW symmetry 
breaking. 

2. IS A JET [19]: the ISAJET 7.69 package contains ISASUSY, which treats the pro- 
duction and decays of supersymmetric particles. The calculations are done within the 
MSSM framework, if the input parameters are provided at the EW scale. The ISAJET 
package provides in addition the possibility to choose the input parameters at the GUT 
scale and then to perform the RG evolution down to the EW scale. 

The generated SUSY spectra have been linked [20] to the DarkSUSY program to calculate 
the relic density. 

3. DarkSUSY [21]: DarkSUSY 4.00 is currently one of the most advanced programs to 
perform DM calculations in a SUSY framework. It computes the relic density of the 
lightest neutralino which is assumed to be the DM particle. The calculations include 
the impacts of resonances, pair production thresholds and coannihilation processes 
[17]. To check the validity of the input parameters the presently known bounds from 
accelerators are included. This package also computes fluxes for a large variety of 
detector types. 

4.1 Differences between SuSpect and ISAJET 

As the implementation of the RG evolution is different in the two programs, small differences 

in the SUSY spectra and couplings arise as shown in FigsIBl and Ur- The mass spectra 

provided by these two programs agree at a level of about 10% for models with mo not much 

larger than my 2 an d not too large values of tan/3. However, the calculated SUSY masses 

can differ by more than a factor of 2 for models in the focus point region or large tan/3 

values [22]. The neutralino masses calculated by the two programs agree quite well, e.g. for 

the Snowmass point SPS4 (m = 400 GeV, my 2 = 300 GeV, tan/9 = 50, p > 0): in FigHthe 

masses of the four neutralinos are shown as a function of A . The LSP mass is independent 

of A and the approximate mass degeneration as well as the A dependence of the two heavy 

neutralinos is clearly visible. However, for the same SPS benchmark point the masses of 

the sbottom squarks agree only on the 10% level between the two programs, as indicated 

in Fig|7| This difference directly affects the relic density calculations because of significant 

sbottom exchange contributions in the neutralino annihilation process. Due to the smaller 

d Models with Aq larger or smaller than displayed are excluded because of problems in the numerical 
application of the RG equations. 
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Figure 6: Neutralino masses as a function of A$ obtained with the SuSpect (left) and the 
ISAJET (right) programs for the Snowmass point SPS4 (m = 400 GeV, my 2 = 300 GeV, 
tanP = 50, fi> 0). 
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Figure 7: Masses of the two sbottom squarks for the Snowmass point SPS4- The masses in 
the left/right plot are computed with SuSpect /IS -A JET. 
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Figure 8: Relic density as a function of A , obtained with the SuSpect (left) and ISAJET 
(right) programs for the Snowmass point SPS4- The blue circles correspond to a calculation 
without any coannihilation processes, while the red points include coannihilations. The green 
shaded region marks the allowed regions by the WMAP data. 



sbottom mass m~ bi calculated by SuSpect, the resulting relic density (FigJHJ) is larger than 
the corresponding ISAJET value. Since the WMAP constraints on the relic density are quite 
strong, the number of models, which survive these constraints are different for SuSpect and 
ISAJET as demonstrated in FigJHl Moreover, the regions in m^, allowed by the WMAP 
data for fixed values of A , are different, too. 

5 Scan of the supersymmetric parameter space 

Assuming CDM to consist exclusively of neutralinos, the cosmological bounds on the neu- 
tralino relic density Q x h 2 imply strong constraints on the mSUGRA parameter space. Since 
the aim of this analysis is to investigate the impact of Aq on the allowed regions in the 
mSUGRA parameter space, we varied Aq in a first step, while keeping the other four pa- 
rameters fixed. For this study we have chosen the Snowmass points [14] as benchmarks. In 
a second step we also varied m and my 2 for different but fixed values of tan/3. Recent 
data for {g — 2) M favour a positive value of the Higgsino mass parameter [23], therefore we 
focused our analysis on \i > 0. A few million mSUGRA models have been generated for 
the remaining four input parameters. We varied m and between and 2 TeV, tan (3 
between 5 and 50 and Aq within ±4 TeV. For the calculations the mass of the top quark has 
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Figure 9: Relic density as a function ofmyi obtained with the SuSpect (top row) and ISA JET 
(bottom row) programs for tanfi = 10, \i > and different values of Aq. The green shaded 
regions mark the WMAP allowed region. The yellow shaded regions show the corresponding 
my2 ranges, which can differ significantly for the two MC programs. 



been set to 178 GeV. 

For all these models the relic density Q x has been computed and required to be within 
0.094 < Q x h 2 < 0.129. In addition, the accelerator constraints adopted from the 2002 limits 
of the Particle Data Group have been taken into account [24]. There are some updates of 
these constraints [25], including the rare decay b — > sj. The accelerator constraints are 
applied to the mass spectra at the EW scale. 

FiglTOl shows that the WMAP strip described in [3] broadens significantly, if Aq is allowed to 
vary. The extension of the allowed regions originates mainly from larger values of mo, which 
are allowed for non- vanishing Aq, while the minimal mo values are obtained for Aq = TeV. 
Large values of rnyi are excluded by the WMAP constraints on the relic density. Since 
slightly larger values of £l x h 2 are derived by SuSpect than by ISAJET, more SuSpect models 
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Figure 10: Allowed models in the m — my 2 plane for tan(3 = 10 and /i > obtained 
with SuSpect (left) / ISA JET (right). The black line corresponds to the WMAP strip for 
A = TeV from Ref. [3]. 



with larger my2 are excluded. The parametrisation of the WMAP strips (black lines in 
FiglTUI) have been determined with other programs [3], thus introducing small differences 
with respect to our analysis. 

In order to avoid colour and/or charge breaking the trilinear scalar couplings at the EW 
scale A t) 6 iT have to be approximately constrained as [26]: 

Al<3(m% u +ml L +ml R ), 

A 2 b <3(m 2 Hd + ml L +ml R ), (5) 
Al<Z(m 2 Hd +m\ L +m 2 fR ). 

The consequences of applying these cuts is depicted in Fig^Jfor tan f3 = 10 and /i > 0. About 
35% of the mSUGRA models, generated with SuSpect, satisfying the WMAP constraints 
are excluded, while for ISAJET the impact of these cuts is even stronger, i.e. about 45% of 
the models are rejected fFig lTTj) . By far the biggest effect originates from the cut on A T due 
to the light stau masses rrif L R . Most of the models (more than 90%) with Aq < — 1 TeV are 
excluded by the cuts. Models with A values larger than 1 TeV are significantly affected, 
too (about 40%). 

To avoid CCB these cuts are necessary but not sufficient, since the vacuum expectation 
values of the squarks, the sleptons and the corresponding Higgs boson were assumed to be 
equal, for simplicity. Moreover these bounds (Eq{HJ) were derived from the tree level scalar 
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Figure 11: Allowed models in the m,Q — my 2 plane for tan/3 = 10 and \i > obtained with 
SuSpect (left) / ISA JET (right) after applying the cuts on A t} b, T to avoid colour and/or charge 
breaking. The black line corresponds to the WMAP strip for Aq = TeV from Ref. [3J. 



potential, while radiative corrections are expected to modify them. 

The scalar potential may contain global CCB minima in addition to the local EW breaking 
minima. As no CCB has been observed, the Universe in its present state may be trapped in 
a local EW breaking minimum. Since this metastable state may have a lifetime longer than 
the age of the Universe due to the small tunnelling probability into the global minimum [27], 
CCB cannot be excluded. Thus, we did not apply these cuts in this analysis. 
In Fig^J all allowed models for different values of the trilinear scalar coupling Aq and tan /3 
are shown in the mo — my 2 plane assuming the Higgsino mass parameter fx to be positive. 
This plot corresponds to FigEl but with tan/3 = 5, 10, 20, 35 and 50 superimposed. In 
contrast to Fig0 where only a few lines survived the WMAP constraints, extended regions 
in the mSUGRA parameter space are allowed, if Aq is varied. For tan/? = 10 mSUGRA 
models with mo ~ 350 GeV for largely negative Aq are within the WMAP constraints, while 
for Aq = TeV the upper bound was about m ~ 200 GeV. The analogous behaviour can be 
observed for larger tan f3 values, where mSUGRA models with mo around 2 TeV are allowed 
for large \Aq\. The masses of the gluinos, of the first two squark generations and of the 
light scalar Higgs boson are nearly independent of Aq. Since the squarks of the first two 
generations are almost degenerated in mass, only m {il is shown in FigE3 as representative 
examples. However, the third generation squark masses 2 , m^ 2 as well as the heavier 
Higgs masses mn t A,H ± depend on the trilinear coupling Aq in general. The future Large 
Hadron Collider (LHC), with an anticipated luminosity of 100 fb _1 , will cover the my 2 
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tang = 5 



500 1000 1500 2000 

m 1/2 (GeV) 



Figure 12: Allowed models in the rao-mj/2 plane for mo ,my2 < 2 TeV, tanf3 between 5 and 
50, fi > and A Q within ±4 TeV from ISAJET. The black lines correspond to the WMAP 
strips as parametrised in [3], with tanf3 = 5, 10, 20, 35 and 50, where tan(3 = 5 belongs to 
the lowest black line and tanf3 = 50 to the highest one. The brown lines indicate the LHC 
discovery reach for an integrated luminosity of 100 fb^ 1 and 300 fb~ x [28], respectively. 

domain up to about 0(TeV), almost independent of mo [28] (FigIT2"J). Thus, most of the 
models generated in this analysis are within the reach of the LHC. 

6 Parametrisation 

As can be inferred from Figs fTH [T5l and EH the allowed regions in the m — my2 plane 
correspond to lines if the trilinear coupling A$ is kept at fixed values. Because of their 
smooth narrow shape, they can be fitted by a polynomial of 2nd order: 

m = a + b ■ mi/2 + c ■ m\j 2 . 

The coefficients of the quadratic terms are small, but they should not be neglected aS C ' TTly2 
can be of the same magnitude as b-my 2 in the allowed domain: the surviving mSUGRA 
models do not lie on a straight line but on a curve fFigs!T^ltdl6|). The parameters we obtained 
by using the MINUIT [29] routines are given in TabsEltolBlfor different values of tan/3 and 
A . The numbers have been calculated from the ISAJET output. 



15 



700 r 
600- 
500- 
400- 
300- 
200- 
100 — 



o r i i 
-1000 




-500 



500 1000 
A (GeV) 



> 



700 r 
600 
500 
400 
300 
200 
100 


























-m H± 






— 








- m H,A 














-m h 






— 


1 










-• 


































































> > , > 





1000 



-500 



500 1000 
A (GeV) 



Figure 13: Gluino, squark and Higgs masses for SPS1A (m = 100 GeV, mi/ 2 = 250 GeV, 
tanf3 = 10, fi > 0) as a function of A obtained with ISAJET. 
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- 990 GeV 
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203 GeV 


Table 2: 


Coefficient 


a, b and 


c of th 


e parametrisation for tan[3 = 


5 for different discrete val- 



ues of Aq between —2TeV and +2.5 TeV. For larger or smaller Aq values too few mSUGRA 
models survive the WMAP constraints to allow a reasonable parametrisation. The last two 
columns contain the domains for m\/2 and tjiq, respectively. 

The gaps on the fitted lines for tan /5 = 50 originate from the chosen step size for m (FigEJ). 
Some plots exhibit extended areas for small m and small (for example Fig^] for 

Aq = TeV). These regions are excluded in the fits, since a proper parametrisation would 
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1000 
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0.000004 
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984 GeV 
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203 GeV 


1500 
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0.016 ± 0.009 


0.000097 ± 


0.000006 
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0.00001 
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Table 3: Same as in Tab\^ but for tan/3 = 10. 
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Figure 14: Allowed regions in the m — mi/ 2 plane, on the left for tan/3 = 5 and on the right 
for tanf3 = 10. In both plots \i is positive. The black lines are the fits given in Tabs\^and\^ 

be much more involved. 

The relic density develops a strong dependence on the common gaugino mass my2 for small 
values near the lower bound coming from the accelerator constraints. Consequently, unex- 
pectedly large m values (for small my 2 ) can lead to Q x h 2 values within the WMAP range. 
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Table 4: Same as in Tab\^ but for 


tanf3 = 


= 20. 
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m 
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0.19 ± 0.02 


0.00004 ± 0.00001 
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799 GeV 
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600 GeV 
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301 ± 2 


0.190 ± 0.005 


0.000051 ± 0.000004 
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-500 


176 ± 1 
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0.000054 ± 0.000003 
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0.000047 ± 0.000002 
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1000 GeV 


151 - 


388 GeV 


500 


138 ± 1 
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0.000131 ± 0.000005 
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0.000122 ± 0.000006 
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0.000103 ± 0.000010 


517- 
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Table 5: Same as in Tab\^ but for 


tanf3 = 


-- 35. 







This behaviour manifests itself e.g. in the isolated points in Fig^](tan — 10) and in FiglTBI 

(tan/5 = 20), for A = -1500 GeV. 

The separation between the lines for different Aq becomes larger with increasing tan (3 and 
increasing \Aq\ fFigsfHI to ITBj). For tan/5 = 5 the corresponding splitting in mo is smaller 
than 90 GeV, while for tan/3 = 50 it can reach about 700 GeV. For negative Aq the shift 
in mo is larger than for positive values, but always to higher mo values so that the minimal 
mo values are obtained for vanishing A . For small tan/3 the dominant contribution to the 
annihilation cross section may come from the t-channel stau exchange (especially in the 
coannihilation region m ^1/2) • The running of the stau masses depends quadratically 
on Aq. For large values of tan (3 the rapid annihilation via a s-channel Higgs boson exchange 
can be dominant. The Higgs masses depend linearly as well as a quadratically on Aq [1]. 
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Figure 15: Allowed regions in the mo — my2 plane for tan/3 = 20 (left) and for tanf3 = 35 
(right). In both plots fi is positive. The black lines are the fits given in Tabs^ and\^ 
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Table 6 




Same 


as in Tab\^ 


I but for 


tan(3 = 


50. 







These dependencies of the sparticles and Higgs bosons masses on A Q fFiglTB^) affect the cross 
sections and therefore also the relic density [30]. They explain the different behaviour for 
positive and negative values of A . 
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Figure 16: Allowed regions in the niQ — mi/2 plane for tanj3 = 50 and \i > 0. The black lines 
are the fits given in Tab\$ 

7 Conclusions 

By including the trilinear scalar coupling A as a free parameter in the range of ±4 TeV 
large regions in the tuq — m-u2 plane of mSUGRA models turn out to be consistent with 
the WMAP data. Fixing this coupling results in narrow lines strongly depending on Aq. 
Large m values (e.g. up to 2 TeV for tan/3 = 50) are allowed for \A \ ^> 0. Using the 
mSUGRA models generated with the ISAJET program we constructed a parametrisation for 
the common scalar mass mo as function of the common gaugino mass for several fixed 
values of A and tan/3. The same qualitative behaviour is found using the SuSpect code. We 
analysed the discrepancies between the two Monte Carlo programs and traced them back to 
slightly different mass spectra and couplings, which result in different values of Q x h 2 . As 
the WMAP cut on the relic density is very tight a small difference may have a big impact 
on the mSUGRA parameter space. 

In addition we started to study the effect of A on the masses of the sparticles. The gluinos 
and the first two generations of the squarks are insensitive to Aq, hence most of the allowed 
mSUGRA models still lie within the reach of the LHC. However, the third generation of 
squarks and the heavier Higgs bosons depend significantly on Aq. 
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